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ABSTRACT 

Radio polarization observations provide essential information on the degree 
of order and orientation of magnetic fields, which themselves play a key role in 
the particle acceleration processes that take place in supernova remnants (SNRs). 
Here we present a radio polarization study of SN 1006, based on combined VLA 
and ATCA observations at 20 cm that resulted in sensitive images with an an- 
gular resolution of 10 arcsec. The fractional polarization in the two bright radio 
and X-ray lobes of the SNR is measured to be 0.17, while in the southeastern 
sector, where the radio and non-thermal X-ray emission are much weaker, the 
polarization fraction reaches a value of 0.6±0.2, close to the theoretical limit of 
0.7. We interpret this result as evidence of a disordered, turbulent magnetic field 
in the lobes, where particle acceleration is believed to be efficient, and a highly 
ordered field in the southeast, where the acceleration efficiency has been shown 
to be very low. Utilizing the frequency coverage of our observations, an average 
rotation measure of ~ 12 rad m~ 2 is determined from the combined data set, 
which is then used to obtain the intrinsic direction of the magnetic field vectors. 
While the orientation of magnetic field vectors across the SNR shell appear ra- 
dial, a large fraction of the magnetic vectors lie parallel to the Galactic Plane. 
Along the highly polarized southeastern rim, the field is aligned tangent to the 
shock, and therefore also nearly parallel to the Galactic Plane. These results 
strongly suggest that the ambient field surrounding SN 1006 is aligned with this 
direction (i.e., from northeast to southwest) and that the bright lobes are due 
to a polar cap geometry. Our study establishes that the most efficient particle 
acceleration and generation of magnetic turbulence in SN 1006 is attained for 
shocks in which the magnetic field direction and shock normal are quasi-parallel, 
while inefficient acceleration and little to no generation of magnetic turbulence 
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obtains for the quasi-perpendicular case. 

Subject headings: acceleration of particles — ISM: individual objects (SN 1006) - 
ISM: supernova remnants — magnetic fields — polarization 
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Introduction 



Shock waves driven by supernova remnants (SNRs) are proposed to accelerat e particles 
to relativistic energies th r ough the di ffu sive shock acceleration (PSA ) mechanism (jKrymskii 



1977 



Axford et al 



1977 



Bel 



1978a 



Blandford fc Ostriker 



19781 ). In this process, also 



known as first-order Fermi acceleration, electrons, protons, and other ions scatter back 
and forth on magnetic irregularities across the shock, gaining energy on each successive 
crossing. The orientation of the ambient magnetic field with respect to the direction of 
the shock velocity and the level of turbulence of the field throughout the shock transi tion 



region play key roles in the DSA process (for a review see 



Malkov fc O'C Drurv 



20011 ). A 



critical issue of the DSA theory concerns the amplification of the ambient magnetic fie 



which is now b elieved to occur synergistically with particle acceleration (jBell fc Lucek 



2001 



Bell 



20041 ) . In this picture, accelerated particles streaming ahead of the shock excite 
magnetohydrodynamic turbulence that can amplify the ambient magnetic field far above 
its in itial seed v alue. The field so produced is expected to be highly turbulent (6B/B > 1) 



.e.g., 



Bell 



20041 ) . Radio polarization observations are therefore crucial to understanding 



such processes, since they can provide information on the degree of order and orientation of 
the magnetic field in the shocked gas of the remnant. Bilateral SNRs are ideal targets for 
these studies — the highly symmetric morphology of such remnants suggests evolution in a 
magnetized environment with large-scale order. 

The remnant of SN 1006, a large (~ 30 arcmin) bilateral SNR located at high Galactic 
latitude, is an ideal laboratory to study the polarization and magnetic field distribution 
while avoiding confusion from other sources in the line of sight. It is widely accep ted that 



SN 1006 is the remnant of a Type la supernova event (IStephenson fc Green 



2002 



p. 174). 



Optical prop er motion measurement s, combined with the shock velocity as inferred from 



Ha spectra (jGhavamian et al. 



20021 ). place this remnant at a distance of 2.18±0.08 kpc 
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( Winkler et al. 


2003 


). At both radio and X-ray wavelengths (e.g., 


Reynolds & Gilmore 


1986; 


Cassam-Chenai et al. 


2008; 


Dver et al. 


2009; 


Rothenflue et al. 


2004 


), the SNR 



appears as a circular she 



contrast, the Ha image (IWinkler fc Long 



1, with two bright lob es perpendicular to the Galactic Plane. In 



19971 ) reveals a system of faint, relatively thin 



filaments defining an almost perfect, complete circular ring, with enhanced emission at the 
NW, between the bright radio and X-ray lobes. Based on XMM-Newton X-ray observations, 



Acero et al. 



(120071 ) estimated a higher ambient density at this location. TeV energy 7— rays, 



originating in the radio and X-ray bright lobes, have recently been detected with H.E.S.S. 



(lAcero et al 



2010|) 



Kundul (jl970j) using 6 and 



The first full polarization map of SN 1006 was obtained by 
11-cm observations performed with the NRAO 43-m telescope at Green Bank. The SNR 
was found to have an average polarization percentage of 10% to 15 % at 5 GHz and a radi al 



Dickel fc Milne! (1 1979 1 



orientation of the magnetic field. A similar result was obtained by 
using higher resolution observations, albeit with a lower fractional polarization (~ 10% 
he h ighest resolution polarization study of SN 1006 was done by 



Reynolds fc Gilmore 



(119931 ). based on VLA observations at 1370 and 1665 MHz. They found peak polarization 
fractions of 30% in the lobes, somewhat higher than previous results, and that the magnetic 
field, although presenting mostly a radial orientation, was predominantly disordered. 



Moreover, they also detected weak polarized flux from the faint SE quadrant. 



ain the bilateral symme t ry of SN 1006 



One of the current theories to exp 


lain t 


Fulbrieht & Reynolds 


1990 




Dver et al. 


2004: 



Petruk et al 



2009 



Schneiter et al 



(e^ 



20101 ) 



interprets the limb-brightened emission of the lobes as a result of a shock expanding into 
a region of the interstellar medium (ISM) whose magnetic field is perpendicular to the 
shock normal. In the case of SN 1006, that means that the magnetic field would be aligned 
SE to NW (the equatorial model). Acceleration rates are higher in perpendicular shocks 
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( lEllison et al.lll995t |Jokipiilll9871 ). as the magnetic field is directly compressed and charged 
particles are accelerated as they traverse magnetic gradients; however, the injection of 
thermal particles is lower due to the high obliquity of fields with the shock normal, requiring 
a mechanism for scattering thermal electrons and ions back into t he shock. However, 



other work (e.g., 



Rothenflug et al. 



2004; 



Cassam-Chenai et al. 



2008) on the variation of the 



non-thermal synchrotron X-ray and radio emission with azimuthal angle around the limb of 
SN 1006 casts doubt on this interpretation. These results suggest that the morphology of 
the emission is best explained if the bright limbs are polar caps, requiring that the ambient 
magnetic field direction extends from SW to NE. Acceleration rates in the bright caps in 
this situation (quasi-parallel shocks) would be comparatively lower than in the former case 
(quasi-perpendicular), but injectio n rates of therma l parti cles would be much higher. The 



simple geometrical argument that 



Rothenflug et al. 



(J200J) applied to the XMM-Newton 



image, favors the equatorial belt model when applied to a radio imag e . To reconcile 



the model with the observations in both spectral ranges, 



Petruk et al. 



(120 111) propose a 



somewhat contrived model of a magnetic field with a gradient. 

When DSA is very efficient, the back reaction of the accelerated particles on the 
structure of the supernova remnant can have significant and observable consequences. One 
of these, potentially active in young ejecta-dominated remna nts, is a narrowing of th e 



distance from the forward shock to the contact discontinuity (Decourchelle et al. 



effect that was first observed in the Tycho SNR ( I Warren et al. 



20051 ) 



2000), an 



Cassam-Chenai et al. 



( 120081 ). using Chandra X-ray observations, discovered this effect in SN 1006 by comparing 
the outermost location of the Ha filaments to the outermost extent of the ejecta, as traced 
by the 0.5-0.8 keV band which is dominated by emission from highly ionized O ions in 
the SN ejecta. This study found that the gap between the forward shock and contact 
discontinuity was largest at the SE rim and that the gap smoothly decreased with azimuthal 
angle as one proceeded toward the bright lobes. The authors concluded that this was direct 
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evidence for a variation in acceleration efficiency as a function of azimuth around the rim of 
SN 1006, likely due to a corresponding azimuthal variation in the orientation between the 
shock normal and the ambient magnetic field direction. 

In this paper we offer insights on these critical issues by carrying out a new high 
resolution polarization study of SN 1006 based on observations performed with the Australia 
Telescope Compact Array (ATCA) and the Very Large Array (VLA). The combination of 
these two data sets allows us to construct images with unprecedented detail, from which 
new information is derived to clarify the role of magnetic fields in the emission of SN 1006. 

2. Observations and Data Reduction 

Full-polarization radio continuum observations of SN 1006 were carried out in 2003 
using the Very Large Array (VLA) of the National Radio Astronomy Observatory and 
the Australia Telescope Compact Array (ATCA). VLA observations were performed on 24 
January 2003 for a four-hour duration in its CnD array configuration (primarily used to 
observe sources at southern declinations). Visibility data were recorded from two 12-MHz 
frequency bands centered at 1370 and 1665 MHz. We used the source 3C 286 as our flux 
calibrator, and 1451—400 as the phase calibrator for SN 1006. ATCA observations were 
performed on three separate occasions: 12 hours in the 6B configuration on 24 January, 
12 hours in the 6A configuration on 3 March, and seven hours in the 750C configuration 
on 12 June. Visibility data were recorded from two 128-MHz frequency bands centered at 
1384 and 1704 MHz, divided into 32 channels each. We used PKS 1934-638 as our flux 
calibrator, and PKS 1458—391 as the phase calibrator. 

Flux, phase and polarization calibration of VLA data was performed using the 
Astronomical Image Processing System (AIPS). The ATCA data was calibrated using the 
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Miriad software package ( ISault et al.lll995l ). We used Miriad to perform the remainder of 
the reduction process. The Stokes visibility data from all observations were split into two 
data sets: 1370 MHz VLA data were combined with 1384 MHz ATCA data in one set, 
and 1665 MHz VLA data were combined with 1704 MHz ATCA data to form the other 
set. We created uniform-weighted 'dirty' images of Stokes I (total intensity), Q and U 
for both frequency data sets using the Miriad task INVERT. Images were mosaicked to 
handle a slight pointing offset between the VLA and ATCA observations. The Miriad task 
PMOSMEM was used to perform a joint maximum-entropy deconvolution of the mosaicked 
dirty images. This program can operate on all Stokes images handling the positive and 
negative emission found in Stokes Q and U maps. Although the resulting resolution of the 
dirty maps was ~ 8" x 6" , we restored them during deconvolution with a circular beam of 
10" to improve their signal- to-noise. The total int e nsity image constructed with these data 



was previously published in ICassam-Chenai et al.l ( 120081 ) . 



We used the Miriad task IMPOL to combine the Stokes /, Q and U to form images of 
linear polarized power (Ip = \/ Q 2 + U 2 ), polarized fraction (lp/J), linear position angle 
(if) = 1/2 tan -1 (?7/Q)), a rotation measure map between frequency data sets, and all of 
their associated error maps. All polarization maps (except total intensity) were blanked 
where Ip was lower than 2a, where a is the noise level for the Stokes Q and U maps. We 
used <tq = <jjj = 0.085 mJy beam -1 . The Ricean bias was removed in first order in the 
same task by computing the polarized power as \JQ 2 + U 2 — a 2 , hence the minimum value 
of the polarized emission is ~ 0.15 mJy beam -1 . In Fig. 1 we present the resulting total 
intensity (1(a)) and polarized intensity (1(b)) maps of SN 1006, where in the latter, pixels 
where the total intensity was lower than 0.25 mJy beam -1 (the outermost contour of SN 
1006 according to Fig. 1(a)) were further blanked. 

The integrated flux density of our total intensity map is ~ 17 Jy, in agreement with 
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Fig. 1. — (a) Total intensity image of SN 1006 at 1.4 GHz. The intensity scale, in mJy 
beam -1 , is displayed on top of the figure. The beam size is 10 arcsec. (b) Polarized intensity 
at 1.4 GHz. The intensity greyscale and resolution are as in (a), although the former with a 
different scale. 



Cassam-Chena'i et al. 



(120081 ). where an independent calibration of the same data was 



presented. This v alue is lower than the combined single-dish and interferometer flux of 



18.68 Jy found by 



Dyer et al 



(120091 ) at the same frequency. We expect a shortage due to 



missing flux from short uv spacings. We also recovered an integrated flux of ~ 2 Jy from our 



linear polarizat ion map after primary 
observations of 



jeam correction, in excellent agreement with the VLA 



Reynolds fc Gilmord (119931 ) . Their images also include central emission in 
excess of that found in our images. Indeed, the integration time, bandwidth and uv coverage 
of differing observations, array configurations, and interferometers can result in images with 
different background noise levels and integrated flux densities, but we must also note here 
that the use of different image processing techniques can too. To verify the integrity of our 
process and resulting images, we processed our VLA and ATCA data separately, and found 
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that the polarization maps were also similar. We are confident the combined reduction 
did not corrupt the resulting polarization maps, and the combination of VLA and ATCA 
visibility data actually improved the overall sensitivity toward small-scale structure, as can 
be appreciated, for example, towards the background source ESO 328-IG 013 (section 3.1). 

As mentioned above, the polarization position angle maps at the two frequencies 
were combined by IMPOL to form a rotation measure (RM) map that could be used to 
derotate the position an gles in order to c erive the i ntrinsic direction of ; he ele ctric field 



Reynolds fc Gilmore 



fll993[ ) show that 



vectors. Earlier work by iDickel &: Milnd (119761 ) and 
the RM toward SN 1006 does not cause position angles to wrap more than 2ir radians, 
so we computed the RM for each pixel using only two frequencies, assuming no 7r radian 
ambiguities. The resulting image is shown in Fig. |2j We found that the distribution of 
RM values across both the NE and SW lobes of the shell is well represented by a Gaussian 
distribution centered at 12 rad m -2 , with a la deviation of 20 rad m -2 . This level of 
deviation is expected, considering that the use of a 2a cutoff in the Stokes Q and U maps 
corresponds to a position angle error of 14 degrees. Inspecting each lobe separately, we 
find that the RM distributions peak at 12 rad m -2 and 14.5 rad m -2 for the NE and SW 
lobes respectively. We found no evidenc e of significant differences of the RM value in each 



lobe, unlike that found in G296. 5+10.0 (lHarvey-Smith et al. 



2010) , another high-latitude 
bi-polar SNR whose RM varies significantly, from —92 in one of the lobes to +44 rad m -2 
in the other one. In G296. 5+10.0, the difference is attributed to a toroidal magnetized wind 
from the progenitor, which is not the case for SN 1006, where the similarity in the RM 
values for both lobes is quite reasonable considering the symmetry of the source. 

Our RMs agree with previous polarization studies within the quoted errors. 
Reynolds fc Gilmore! ( 119931 ) found an RM of 28 ± 16 rad m" 2 on the NE lo be, and 16 ± 14 



rad m 2 on the SW lobe. A single-dish polarization study by 



Milne! (119711 ) at 2.7 and 5.0 



■~ -43°00'[ 



Right Ascension (J2000) 



Fig. 2. — Distribution of the rotation measure on SN 1006, at 10 arcsecs resolution. The 
scale bar, shown at the right, is given in radians per square meters. 

GHz yielded RM values of 21 ± 5 rad m~ 2 in the NE, a nd 9 ±5 rad m -2 in the SW, although 



a re-analysis of similar data by 



Dickel fc Milnd (119761 ) had an error of 23 rad m 2 . While 



the RM appears different on average between the two lobes, they are not significantly so. 
Using the assumption that our average RM value of 12 rad m~ 2 was uniformly distributed 
over the remnant, the lowest frequency polarization position angle map was derotated to the 
position angle at zero frequency, and further rotated 90° to display the orientation of the 
magnetic field. The resulting magnetic vector map is shown in Fig. [31 Pixels were blanked 
where the flux density fell below 10 mJy beam -1 in a total intensity image convolved to 60 
arcsec. The magneti c vectors show an oyerall ra dial direction, and appear very much the 



same as reported by 



Reynolds k Gilmore! fll993h . 



3. Results 
3.1. Fractional polarization 



Figure H] shows the distribution of the fractional polarization p, computed as described 
in the previous Section. The p map was compared to the error map a p , and those 
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pixels for which the ratio p/a p < 2 were blanked. In this image, unpolarized emission 
is represen ted in red. The average fra ctional polarization in the SE and NW lobes is 



17 ± 7% 



o . 



Reynolds &: Gilmord (119931 ) note that since their images do not include single 
dish observations, a smooth component of the total intensity is missing, so that if all the 
polarized flux is assumed to be detected, the fractional polarization measured is only an 
upper limit. They discuss this problem and conclude that the ratio maps must show values 
very close to the true ones, especially near the rim. Since the total flux we measure is much 
closer to the real one, the fractional polarization should be accordingly closer to its true 
value. 



Reynolds fc Gilmord (119931 ) called attention to a faint arc on the SE rim visible in 
the Stoke parameter Q map that closely follows the circular rim of the shell. The authors 
argued that this feature is real, however their sensitivity did not allow them to provide 
additional information about its morphology or polarization properties. We have examined 
our Q image and found that the highest emission occurs at an arc that extends from 113° to 
180° at a radius of 880 arcsec centered at RA(2000)= 15 h 2 m 53"36, Dec.(2000) = -41° 56'25". 
Between 113° and 140°, the emission at the arc is several times above the noise, while the 
rest of the feature is weaker. In Fig. [5] we show the average emission in the Q parameter on 
4-arcsec wide annuli from the center over the sector defined by the full arc. The azimuthal 
average clearly enhances the emission at the arc, confirming that this feature is real. The 
same result is observed when only the weakest part of the arc is considered (dashed line in 
Fig. [5]). In Fig. HI this feature appears to be highly polarized, with p = 0.6 ± 0.2, where the 
uncertainty is the RMS of the observed polarization fractions rather than the error on the 
mean. High polarization is also observed at the NW rim. These results will be discussed in 
the next section. 

Finally, these observations allowed us to measure the polarization of a bright 



13 



extragalactic source that lies on the remnant's shell at a position angle of about 9 0° . 



source is cataloged in NED as ESO 328-IG 013, an interacting galaxy (ILauberts 



19821 ) 



'his 



it is also coincident with a 2MASS galaxy (2MASX J15040364 -4155511) at a redsh i 



Reynolds k Gilmorel jl993 ) 



t of 



z = 0.039 and is also likely the radio source PMN J1503-4145. 
noted that the fractional polarization was too high to be interpreted as the core of a radio 
galaxy or quasar. In fact, they found a fractional polarization of about 11% because they 
could not resolve the core of this elongated source. Our data reveal two polarized jets 
emerging from a clearly unpolarized core (Fig. [6]). The jets extend for at least 2' (~100 kpc) 
and have typical polarizations of 18%±5% on the shell of SN 1006 and 65%±15% beyond 
the rim, and the core has p = 3%±0.8%, m ore in line with the value expected for the core of 



an extragalactic source (e.g., 



Saikia 



19991 ). The high polarization of that portion of the jet 
that lies beyond the rim is probably close to its true intrinsic value; however, over most of 
its extent the jet is superposed on the SNR emission which makes the fractional polarization 
appear lower than its true value. To support this argument, we note the remarkable 
difference between the flux density of the galaxy just to the East of the edge of SN 1006 
(~ 0.65 mJy beam -1 ) and immediately to the West (~ 2 mJy beam -1 ), while the polarized 
flux for the same regions is about 0.46 and 0.52 mJy beam -1 respectively. Moreover, the 
flux of SN 1006 just to the North and South of the galaxy in the inner border of the shell 
varies between 1.3 and 1.5 mJy beam -1 , implying that there is a major contribution from 
the SNR to the flux measured on the galaxy. The contribution from the SNR decreases 
slowly along ~ 1 arcmin inwards from the SNR edge, and much faster as it gets closer to 
the core of the galaxy. At the core, the flux density of the galaxy is 20 times higher than 
the SNR. A detailed study of the polarization measurements on ESO 328-IG 013 will be the 
subject of a separate paper. 
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3.2. Magnetic Field Orientation 



Figure [3] shows the distribution of magnetic vectors afte r combining the ATCA and 



Reynolds fc Gilmore! (jl993[h 



VLA polarization observations. A comparison with Fig. 8 in 
where only VLA observations were used, shows a close agreement in those regions where 
data are available, although they detected higher levels of polarized emission in the center of 
the remnant. Based on previous studies, in which polarized emission could only be detected 
on the brightest lobes, it was suggest ed that the magnetic fi elds w ere predominantly radial. 



There were, nevertheless, hints from 



Reynolds &: Gilmord (119931 ) of a non-radial field 



orientation along the faint SE rim based on two, barely significant, magnetic field vectors. 
Our data allow us to gain insight on the details of the magnetic field distribution at the 
outer rim and at the lowest intensity regions in the SE and NW. 

In order to determine whether the orientation of the magnetic field is radial along the 
whole SNR, we computed the angular difference between each magnetic vector and the 
radial direction, a polar-referenced angle , by adopting a polar coord inate system centered 



Reynolds fc Gilmord (119861 ) and 



Winkler fc Long 



at the SNR's nominal center, located by 
(119971 ) at RA(2000)= 15 h 2 m 51!7, Dec.(2000) = -41° 56'33". The polar-referenced angle 



varies between 90° and —90°; it is negative for a clockwise shift with respect to the radial 
direction and positive for a counter-clockwise shift. Since the two extreme angular values 
represent the same orientation (tangential), we used a color scheme that allows us to 



visualize both with the same color. In the plot shown in Fig. 7(a) yellow and green colors 



correspond to radial magnetic field orientations, while blue colors indicate tangential ones. 
It is clear that the overwhelmingly radial direction in the bright lobes is not preserved in 
the SE and NW, where the orientation is predominantly tangential. Even in the lobes, 
there are some departures from the radial direction, mainly at the rim. 



Furthermore, we selected only those pixels with radial field orientations within ±14°, 
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which is the error in the position angle (Section 2), and plotted them in green in Fig. 



7(b) We also selected vectors parallel to the Galactic Plane, which lies at a fixed angle 
of 60° (counterclockwise from North) as projected in this image, and overplotted them in 
red. In both cases, the scale is such that pixels in the limits of the tolerance interval are 
displayed with the lowest intensity. Finally, pixels not belonging to any of these two groups 
are overplotted in blue. In this case, the faintest pixels correspond to angles closest to the 



Galactic Plane direction. Along this direction (yellow line in Fig. 7(b)), pixels appear in 



yellow, clearly indicating that the magnetic vectors there fulfill both criteria. 



We examined the distribution of polar-referenced magnetic field direction (Fig. 7(a) ) 
solely using information from pixels within the NE, SW and SE regions of SN 1006. We 
extracted the information from the bright lobes using pixels within two radial wedges with 
azimuth ranges from 20° to 100° (NE) and 170° to 250° (SW), and radial ranges from 540 to 
860 arcsecs from the center. We found that the distributions peak at 9° in the NE lobe and 
-4° in the SW lobe, with standard deviations, er, of 20° and 25°, respectively. We extracted 
information from pixels in the SE using a radial wedge with an azimuth range from 113° to 
147° , and a radial range from 865 to 920 arcsec from center. We shifted the difference angle 
range from [—90°, 90°] to [0°, 180°], a more appropriate range for identifying magnetic fields 
that are tangential (purely tangent would then be 90°). The distribution in the SE peaks at 
~ 112°, with a standard deviation of 18°, clearly consistent with a tangential orientation. 

A surprising result is obtained if the outer radius of the SW lobe sector is extended 



so as to include the outer rim (see Fig. [8]). In Fig. 8(a), the logarithmic distribution of 
polar-referenced field angles, including only the inner sector (radius of 860 arcsec) is plotted 
in red, while the green line displays the distribution at a radius of 920 arcsecs. While the 
distribution is the same for polar referenced angles between approximately —50° and 50°, 
there is a clear excess of more oblique vectors when the region with radii between 860 and 
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920 arcsecs is included. To further highlight this trend, in Fig. 8(b) we plot the polar 
referenced angle distribution from only the radial range between 860 and 920 arcsec (with 
the distribution shifted to the range [0°, 180°] as before). A simple Gaussian fit (albeit 
crude given the evident non-Gaussian nature of the distribution) gives a peak at ~ 88° 
with a deviation of 30°. Although more sensitive polarization observations are desirable to 
confirm this result, it is worth noting that indications of a ta ngential magnetic fie ld at the 



outer edge have also been reported in the Cassiopeia A SNR (jGotthelf et al 



200lh . 



Finally, we examined the distribution of magnetic field orientation with respect to the 
direction of the Galactic Plane. The distribution, shown in Fig. [91 can be reasonably well 
fitted by the sum of a constant term and two Gaussians (red, dashed lines): a broad one 
centered at 45° and a narrow one centered at 58°. The result is represented by the green 
line. The vertical solid line overlaid at 60° corresponds to the position angle of the Galactic 
Plane in this direction. In Fig. [TOj we consider only the two bright lobes, which provide a 
similar result: 47% of the magnetic vectors are clustered around 58°, while the rest fall in a 
broad component centered at ~ 40°. 



4. Discussion 



Milnd (119871 ) reviewed observations of polarized radio emission for 27 SNRs and by 



examining their projected magnetic fields, concluded that young remnants have radial 



magnetic fields wh ile well-defined tangential fields are seen on 



observations (e.g., 



Reynoso et al. 



1997 



Landecker et al. 



1999 



y in older SNRs. S ubsequent 



Wood et al. 



20081) have 



supported this picture. In the case of SN 1006, the magne tic field orientation could only be 



measured on the bright lobes ([Reynolds fc Gilmore 



1993 



and references therein) and, as 



expected for a young SNR, was found to be radial. These studies left open the question of 
whether fainter regions of the shell had a radial field as well. 
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If the projected magnetic field were radial all around the shell, then the measured 
position angles with respect to North should be uniformly distributed between —90° and 
90°. However, considering that most of the information on the magnetic field in SN 1006 
comes approximately from the first and third quadrants (the NE and SW bright lobes), 
there will be a bias in the distribution favoring position angles between 0° and 90°. While 
the broad Gaussian centered at 45° in Fig. [9] is not conclusive support for this model, the 
narrow component near 60° is clearly at odds with it. When we reproduce the histogram 
for azimuthal sectors of 10° width over the lobes, the fixed angle component near 60° is 
found in most, but not all, of the sectors. 



Examination of Fig. 7(b) clarifies the situation, showing clearly that the pattern of 
magnetic field vector orientations is not uniform. Patches with a radial orientation (colored 
green) are found all along the northern lobe and partially in the southern lobe. On the 
other hand, regions with vectors parallel to the Galactic Plane (colored red) are found 
over practically all of the shell. (The pixels colored yellow are located where the local 
radial direction happens to match the fixed field direction.) The patches with radial and 
fixed field directions are interspersed throughout the bright lobes. There is no correlation 
between the magnetic field orientation and the fractional polarization in these patches. At 
the outermost edge of the remnant's shell, several azimuthal sectors display orientations 
(coded in blue) which are neither radial nor parallel to the Galactic Plane, possibly pointing 
to the addition of a radial and a fixed Galactic component. 

We have established that an important fraction of the magnetic field vectors in SN 
1006 are parallel to a fixed angle near 60°. This is particularly clear in the SE sector, where 
the magnetic field vector distribution is strongly dominated by a fixed field component at 
an angle of ~ 54° and a broadening width of a ~ 10°. We draw the conclusion that the 
fixed angle component indicates the direction of the ambient magnetic field projected onto 
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the plane of the sky, which in this case runs nearly parallel to the Galactic Plane. This 
allows us to answer the long-lasting question concerning which of the two competing models 
explains the bright lobes, i.e., whether they arise from a limb-brightened equatorial belt or 
two polar caps. Our study finds that the two bright lobes in SN 1006 are actually polar 
caps. 

The linear polarization fraction of an optically thin source emitting synchrotron 
radiation in a uniform field is limited by the following relation: p < (a + l)/(« + 5/3). The 
mean spectral index for SN 1006 is a = 0.6, so that the expected maximum value for p is 
70.5%. Such high values are not observed for SNRs. In general, p is me asured to be 10% 



15 %, with higher values (35% - 50%) in Just a few exceptional cases (ILandecker et al. 



1999 



Dickel et al. 



2000 



Reynoso fc Green 



20071 ). The measured fractional polarization 

arization studies 



Burn 



(]1966|), this 



of the SW and NE lobes is about 17%, in agreement with previous po 
towards SN 1006 but far from the theoretical maximum. According to 
means that only 24% of the magnetic energy is in an ordered component. Therefore, even 
though the magnetic fields appear ordered at the level of our image resolution, the low 
fractional polarization indicates more generally that the fields on smaller scales are highly 
disorganized in these regions. 

In contrast, the SE region of faint radio and non-existent X-ray synchrotron emission 
shows a high fractional polarization compatible with the theoretical expected value. This 
requires the emiss ion region to host a highly o rdered magnetic field, likely just the swept-up 
ambient field (see 



Cassam-Chenai et al 



20081 ). Figs. [3]and 7(a) show that the magnetic 



field is mostly tangential at this faint region. This orientation was derived assuming a 
uniform RM for the whole SNR. Were the RM half the average value in this region, the 
vectors would rotate counter-clockwise by ~ 10°, but would still be close to the external 
magnetic field direction estimated from the bright lobes. 
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We can therefore link efficient particle acceleration and disordered, amplified magnetic 
fields with quasi-parallel shocks in the bright SW and NE lobes of SN 1006. Likewise, 
we link inefficient particle acceleration and highly-ordered swept-up magnetic fields with 
quasi-perpendicular shocks in the faint SE rim. This scenario is also in agreement 



with the recent det ection of very high energy 7-rays at the lobes (jAcero et al. 



Petruk et al. 



20101). 



(120121 ) consider the coincidence between the radio, X-rays and 7-ray emission 
as possible evidence for a leptonic origin of the 7-rays, and set upper limits to the magnetic 
field strength of a f ew hundred micro- gauss, which are less constraining than previous 

and references therein). 



determinations (see 



Petruk et al 



2012 



The radial magnetic fields usually observed in young SNRs have been ascribed to 



amplification along the Rayleigh- Taylor finger s developed at ; ; 



shocked ambient medium and the ejecta (e.g., 



Jun fc Norman 



re interface between the 



19961 ) . However, we have 



found that the radial distribution is only partially observed in the bright lobes of SN 1006. 



Dolginov et al 



(120041 ) solved the nonlinear problem of the generation of a magnetic field 
by particles accelerated in a MHD shock front. Their model is based on the assumption 
of uniform stationary particle injection, which is probably suitable for SN 1006 due to the 
low density and homogeneity of the ambient medium, far as it is from the Galactic Plane. 
They found that when the ambient magnetic field has a component perpendicular to the 
shock front, the parallel component is destroyed by the secondary field generated by the 
accelerated particles, resulting in a radial field distribution. Our results for the bright lobes 
suggest that this rearrangement is not complete and a fraction of the original field survives 
as inferred from the peak near 60° in the distribution of magnetic field dir ections. On the 



other hand, if the shock front were parallel to the ambient magnetic field, 



Dolginov et al. 



( 120041 ) find a second possible solution in which the field is not suppressed but, instead, is 
amplified and can be two orders of magnitude stronger than the mean value measured in the 
cold and warm phases of the ISM. We offer the following comments, based on the properties 
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of the SE sector of SN 1006 , wher e the observed field is nearly parallel to the shock front 
First 



Cassam-Chenai et al 



( 120081 ) showed that the faint radio surface brightness there 
could be accounted for by a small injection rate of particles into the acceleration process 
(corresponding to < 10% of the shock energy going into accelerated particles) and shock 
compression by a factor of 4 of a typical 3 /zG ambient ISM magnetic field. In other words 
no magnetic field amplification appears to be required in the SE sector. Second, the high 
degree of radio polarization we find in the current study requires that the post-shock 
magnetic field be highly uniform. The interpretation we favor is that particle acceleration 
and the generation of turbulent magnetic field are inefficient at the SE rim of SN 1006. 

We take the fixed angle component observed almost all along the SNR shell as a strong 
indication of the ambient magnetic field direction. The prevalence of this direction at the 



highly polarized SE arc, which can be explained by the second solution of 



Dolginov et al. 



( 120041 ). further supports this hypothesis. The patchy distribution of regions with either 
fixed angle or radial magnetic vectors, which is a surprising result of our study, remains to 
be understood. 



5. Conclusions 

We have presented a new radio polarization study of SN 1006 that has unveiled new 
features related to the polarization fraction and magnetic field direction of both the SNR 
and the surrounding ISM. These findings have implications for DSA theory. We found that 
the brightest radio, X-ray, and TeV 7-ray features have the lowest polarization fractions, 
indicating the presence of a disordered, turbulent magnetic field precisely in the regions 
where particle acceleration is most efficient. Moreover, the SE region of the shell, with 
very faint synchrotron emission, shows a fractional polarization close to the theoretical 
maximum value for an optically thin source, indicating a highly ordered magnetic field at a 
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20081 ). While 



location where particle acceleration is very inefficient (jCassam-Chena'i et al. 
most of the SNR shell appears to have a radial magnetic field distribution, we have found 
that a significant portion of the distribution share a single angle of ~ 60°, which in turn 
coincides with the Galactic Plane direction. In the SE sector the distribution of magnetic 
field directions is consistent with a single component at an angle of ~ 55°. Overall, the 
evidence points to an ambient magnetic field roughly parallel to the Galactic Plane, of 
which the SNR retains some knowledge even after the passage of the shock front. The 
portions of SN 1006 where the shocks are quasi-parallel are more efficient at producing 
particle acceleration and generating magnetic turbulence; elsewhere where the shocks are 
quasi-perpendicular particle acceleration is inefficient and there is little turbulent magnetic 
field generation. Finally, given these results, considerations of geometrical symmetry clearly 
establish that, rather than being the limb-brightened edges of an equatorial belt, the bright 
synchrotron emitting lobes in SN 1006 are polar caps. 
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Fig. 3. — Distribution of magnetic field vectors on SN 1006 at 1.4 GHz corrected for Faraday 
rotation (assuming uniform RM=12 rad m -2 ), at 10 arcsecs resolution. Total intensity 
contours at 10, 20 and 50 mJy beam -1 , are superposed, where the beam was convolved to 
60 arcsecs. For the vectors, a length of 30 arcsecs represents 0.25 mJy beam -1 of polarized 
flux. 
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Fig. 4. — Fractional polarization p of SN 1006 at 1.4 GHz. The resolution is 10 arcsecs. The 
color scale is shown at the right. Only pixels where p was at least twice its error were kept. 
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Fig. 5. — Azimuthal average of the Stokes parameter Q computed between 113° and 180° 
over annuli of 4 arcsec width. The dashed curve shows the same but averaged between 140° 
and 180°, where the feature at 880 arcsec is weaker. 
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Right Ascension (J2Q00) 



Fig. 6. — Fractional polarization of the extragalactic source ESO 328-IG 013, located on 
the NE lobe of SN 1006. The resolution is ~ 8" x 6". A few radio continuum contours are 
included to show the outer rim of the SNR. 
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(b) Radial and fixed angle distributions 



Fig. 7. — (a) Magnetic field orientation with respect to polar angle (polar-referenced angle). 
The center of the polar coordinate system used to define the polar angle (local radial direc- 
tion) is marked by a yellow cross at the center of SN 1006. The color scheme of the legend 
is cyclic; blue represents both 90° and —90°. A positive polar-referenced angle indicates a 
counter-clockwise angular difference between magnetic vectors displayed in Fig. [3] and the 
polar angle, (b) Magnetic field orientation with respect to the Galactic Plane and polar 
angle. Red pixels are for vectors at a fixed angle of 60° (the direction of the Galactic Plane), 
while preen indicates vectors that are locallv radial. In both cases, a, tolerance of +14° is 
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(a) (b) 

Fig. 8. — (a) Logarithmic distribution of the direction of magnetic vectors in the SW sector 
based on a polar coordinate reference system, extracted from pixels within a radial wedge 
defined between 170° and 250° in azimuth, and from 540 to 860 arcsec in radius (red) or to 
920 arcsec in radius (green), (b) Difference polar- referenced angle distribution (see text). 
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Fig. 9. — Distribution of the position angle of the magnetic field measured with respect to a 
frame fixed to the celestial coordinates (green curve). The best fit is represented by a black 
solid line, which is the result of the addition of two Gaussian distributions (red, dashed lines) 
plus a constant. The vertical solid line corresponds to the Galactic Plane tilt in the direction 
of SN 1006. 
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B position angle (degrees) 



Fig. 10. — Same as Fig. |9]but considering only the two bright lobes. 
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